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ABSTRACT: In this paper we discuss more wide applicability of the Chan- 
drasekhar model of coupled oscillators in the optical rotatory dispersion of the 
crystals. We solve the problem using the Chandrasekhar model of two coupled 
oscillators in the case when we include all couplings between adjacent oscillators 
on the helix that is given by the crystal structure. Further we discuss the re- 
sults of coupled oscillators models in the case of the including of the couplings 
between even and odd oscillators on the helix. The ORD results obtained after 
the approximations of the oscillator strengths verify that these couplings have 
the important influence on the crystal ORD. 

1 Introduction 

The phenomenon of the optical activity (OA) has two aspects which arise from 
the interaction of the radiation with matter - dispersive and absorptive. These 
aspects are the optical rotatory dispersion (ORD) and the circular dichroism 
(CD). It is well known that ORD is the dependency of the rotation of the linear 
polarized light per unit length on the frequency u or on the wave length A. The 
CD is the ellipticity per unit length of the wave getting off the crystals. In the 
difference from the ORD, which is nonzeroth in a wide frequency region, the CD 
is nonzeroth only in a very narrow frequency region in the absorption region. 

The important group of optically active crystals are the crystals with screw 
axis of symmetry belonging to the space groups of symmetry D% and D®. The 
typical representatives of these crystals are a-quartz, cinnabar, tellurium, se- 
len, camphor, benzil. The optical activity of these crystals due to asymmetrical 
originating of the crystal structure because the molecules or atoms forming the 
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crystals are not optically active. Camphor is the exception because his molecules 
are optically active so that his optical activity has two sides - molecular and 
crystalline. 

In the past the OA of crystals was studied by more authors and their works 
are based on the different theories. We can introduce the theory of excitons 
[|], ||, |3], f| , the theory of coupled oscillators |5], |, |7|] or Lagrangian formalism || . 
But it is known that the theory of coupled oscillators gives the results which are 
very good applicable in the fitting of the experimental data of OA. The model 
of coupled oscillators was for the first time used by Chandrasekhar [g| which 
has applicated the Kuhn model of two coupled oscillators which represent 
the smallest unit of the optical active crystal. He has used the following model: 
The first oscillator lies in the plane z = and its position is given by direction 
cosines a, /3, 7. The second oscillator lies in the plane z = d and its direction of 
vibration is turned by the angle 9 around the z axis which is parallel to the crystal 
axis c. Both oscillators lie on the helix which is given by the crystal structure. 
The z axis (crystal axis o) has the direction of the propagation of the linear 
polarized electromagnetic wave. Both coupled oscillators forming one compound 
oscillator are identical. We denote the coupling constant between oscillators as 
Qi. Chandrasekhar further has assumed that the number of compound oscillators 
in the volume unit is N/2, where N is the number of single oscillators in the 
volume unit. The influence of the interactions between oscillators belonging to 
the other helices is neglected because of another type of couplings. 

Becouse of the interaction between single oscillators in the compound oscil- 
lator the natural frequency u is split into two frequencies uj\ and u 2 of normal 
modes of vibrations. For these normal modes Chandrasekhar has solved the 
dispersion theory of refractive indices for the propagation of the left and the 
right circularly polarized wave into which the linear polarized wave is split in 
the optically active medium. It is well known that the medium containing cou- 
pled oscillators is optically active and therefore it must be characterized by the 
different refractive indices n\ and n r for the left and the right circularly polar- 
ized wave. Chandrasekhar has solved the dispersion theory semiclasically, he 
has introduced into results the oscillator strengths but he has assumed that the 
oscillator strengths of the normal modes of vibrations are the same. 

Chandrasekhar has solved his model of optical activity for the a-quartz. Be- 
cause the CD is not measured for the a-quartz and therefore nor the ORD in the 
absorption region Chandrasekhar has assumed the oscillators as undamped. 

Later V. Vysin JF] has removed Chandrasekhar's simplifications in the ap- 
proach of the oscillator strengths of the normal modes of vibrations and he has 
obtained for the ORD which is denoted by p(uj) the formula 



where f qi and f q2 are the oscillator strengths of the normal modes of vibrations, 



TiNde 2 uj 2 
mc 2 



(a 2 + fi 2 ) sinO 
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m is the mass of the single oscillator, e his electric charge and c is the velocity of 



the light. The frequencies of the normal modes are cuf = to 2 + Q 



Qi 



Because the splitting of the frequencies of the normal modes is very small we can 
([l]) rewrite in the form containing the natural frequency of single oscillators 
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We see that the ORD based on the two coupled oscillator model leads in 
general to the two member formula. The first term on the right side of eq. @ 
is known as Drude's term and the second term is known as Chandrasekhar's 
term. This is problem because we don't know the oscillator strengths for the 
real crystals. We can solve the oscillator strengths only in approximations. For 
example in the approximation of the linear harmonic oscillator f qi = f q2 = fo (the 
Chandrasekhar approximation) the formula (0) leads only to the Chandrasekhar 
formula 



p(u) 
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(a 2 + /3 2 ) sin^- 
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On the other hand we can solve the oscillator strengths in the Heitler - London 
approximation [|K|] which gives the result 



fqi 
UJx 

and then ORD is given by the formula 
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which was derived 



We note that the formula (|5]) is the Agranovich type [l 
by the exciton's theory. 

The formula (0) or its special forms (|3]) and (0) was used in many works 
to describe the experimental data of ORD of the crystals, for example of the 
tellur [0, [12|, PI , of the a-quartz [L5| or of the crystal Bii 2 Ge02o @, |17 
and without exception with good results. On the other hand we think that the 
Chandrasekhar model still contains some simplifications which we want to solve 
in this paper. 

The Chandrasekhar model assumes that N/2 coupled oscillators, where iV is 
the number of single oscillators, exist in the volume unit. But it means that we 
have N/2 isolated compound oscillators in the volume unit. This idea desribes 
for example the system of randomly oriented isolated non-interacting molecules. 
But in the crystals all adjacent molecules or atoms interact on the helices. It 
means that the second oscillator in the any compound oscillator is at the same 
time the first oscillator in the next compound oscillator. It is the first question 
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how the results of the Chandrasekhar model will change in the case when we 
include all couplings between all adjacent oscillators into this model. 

The second question follows from the fact that the two coupled oscillators 
model neglect all other couplings between oscillators on one helix. It is possible 
in the case when all oscillators lie on the line because the value of the coupling 
constant decreases with third power of the distance. But it is not true in the real 
crystals where the diameter and the step of the helix are comparable. In these 
cases we must include also the couplings between even and odd oscillators on the 
helices. Further in the crystals with the space groups of symmetry D\ and -Df the 
vibration directions of the adjacent oscillators contain the angle 9 = 120deg. The 
vibration directions of the even or odd oscillators contain the angle 9 = 240deg. 
From the formulae (0), (Q) or (^) we see that the sense of ORD depends on the 
value of sin#. It means that the couplings between even and odd oscillators on 
the helices have the opposite effect on the ORD in compare with the couplings 
between adjacent oscillators. We note that in the crystals with space groups of 
symmetry D% and the couplings between the first and the fourth oscillators 
which are randomly oriented don't have the influence on the OA and then the 
other couplings are really neglectable. 

2 The influence of the couplings between all ad- 
jacent oscillators 

We can discuss this problem by means of the results of our previous paper 
[ y| . In this paper we have solved the application of the ORD Chandrasekhar 
model of three coupled oscillators. We have added to the Chandrasekhar model 
of two coupled oscillators the third oscillator which lies on the helix in the plane 
z = —d and its vibration direction is turned by the angle — 9 around z axis with 
respect to the oscillator in the plane z = 0. We include also as Chandrasekhar 
only the couplings between adjacent oscillators. The number of these compound 
oscillators in the volume unit is N/3. 

The solving method is also the same as Chandrasekhar's with the exception 
that the natural frequency ujq of each oscillators is split into three frequencies uo\ , 
u>2 and CU3 and also the ORD result contains three terms regards to (|IJ) that is 



where uj\ = uj^ + y/2Q\, lo\ = cUq and uj\ = Uq — V2Q\. We can rewrite this result 
with the natural frequency of oscillators u)q but this result is involved too. We 
express this result only in the approximations of the oscillators strengths. In the 
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Chandrasekhar linear harmonic oscillator approximation f qi 

get 



f<n = fga = fo we 



8irNde 2 Qif uj' 2 
3mc 2 



(a 2 + /3 2 ) sin#- 



2Qi cos 9 



(7) 



We see that the second term in the square brackets is very small with respect 
to the first one. It contains the coupling constant Q\ that is very small and the 
denominator of this term is in the frequency region far from absorption much 
greater than the denominator of the first term. From this reason the second term 
can be neglected. 

In the Heitler - London approximation that is given by the extended relation 
(W) we have 



A7TNde 2 Q 1 f {g 2 + (3 2 )smd 
3mc 2 



u z (u$ + a/) 



2Qi cos 9 



or 



(8) 



and we can neglect again the second term in the square brackets. But then the 
results of two and three coupled oscillators are very similar, they differ only in the 
multiplicative constant 4/3. We see the sense of this constant if we compare the 
results for one compound oscillator. We must divide the results of two coupled 
oscillators model (|3]) and (|5]) by N/2 (the number of compound oscillators in the 
volume unit) and the results of the three coupled oscillators model (0) and (||) 
by N/3. By the comparing of corresponding results we see that one triad of the 
coupled oscillators has the twofold effect on the ORD than one couple. From this 
it follows that the effects of couplings between adjacent oscillators are aditive. If 
we want to include all couplings in the two coupled oscillators model we must take 
as the number of coupled oscillators N (in the case N — ► oo) because the number 
of couplings is the same as the number of single oscillators. In the case of three 
coupled oscillators model we take N/2 as the number of compound oscillators in 
the volume unit etc. Then both models of two and three coupled oscillators give 
the same results in the case of the including of the couplings between adjacent 
oscillators only. The mistake of this conclusion is given by the neglecting of the 
second terms in the square brackets in ([?]) and (§) but we can assume that in the 
cases of small couplings in the crystals the mistake is neglectable. 

From our conclusions in this section it follows that if we use our results in the 
approximation of experimental data of ORD of the crystals we get different values 
of unknown parameters in these formulae. For example the ORD of a-quartz 
was approximated by Chandrasekhar's and Vysin's results of the two coupled 
oscillators model. The unknown parameters were the oscillator strengths of the 
normal modes of vibrations. The values of the oscillator strengths were in the 
results that are assumed as true smaller as the number of valence shell electrons 
1{J. In the other results the values of the oscillator strengths were greater. But 



after our correction of the results of the two coupled oscillators model we get for 
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all oscillator strengths a half of values and it would be difficult to say what results 
of the approximations of the experimental data are true and what the values of 
all oscillator strengths will be too small. It is the next reason for the assertion 
that we must include also the couplings between even and odd oscillators which 
due to its opposite effect on the ORD return for example in the case of a-quartz 
the results of the approximations of oscillator strengths to the values that they 
have in ||19||. 



3 The influence of the couplings between even 
and odd oscillators 



The smallest system in which we can include the coupling between even or 
odd oscillators in the helix is the system of three coupled oscillators if we include 
also the coupling between the first and the third oscillator. The scheme of this 
system is the same as in the previous section. We characterize the coupling 
between the first and the third oscillator by the coupling constant Q 2 - Regarding 
to the previous section we use now the value N/2 as the number of compound 
oscillators where N is the number of single oscillators in the volume unit. It 
means that the last single oscillator in any compound oscillator is the first single 
oscillator in the next compound oscillator on the helix. 

The solving ORD of this system by the Chandrasekhar semiclassical model is 
described in 



20]. We get for the ORD the three member formula 
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JQl f q2 COS 6 



UJ 



oj z 
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UJ 2 — id 2 
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-(2 cos 9+ A3) f 
UJ 2 — UJ 2 



(9) 



Q2 + y/8Qf+Q 



LOT, 



Q 2 and uj 2 = uj 2 + 9^M+9l and Al , 



where uj\ = uJq + 

A 3 are A\ = \xh^ — > ^3 = 2qT^ — ■ ^ e see ^ na ^ "J ^ ne expression 

of this result regarding to the natural frequency we obtain complicated result 



and its form is also described in |20|| . Let us again be interested in the form of 
the result after the approximations of the oscillator strengths. If we neglect as in 

we get in the Chandrasekhar 



-oj 



2\3 



previous section all the terms containing 
approximation the result 

AnNde 2 f (a 2 + f3 2 ) sinfl ■ (Q 1 + 2Q 2 cos 



and in the Heitler - London approximation 

2nNde 2 U (a 2 + (3 2 ) sin 6-(Q 1 + 2Q 2 cos 1 



p(u) 
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(u 2 + u 2 ) 



UJq [Uq 
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If we compare these results with the results (|7|) and @ after their discussion, 
it means after neglecting the terms with \/{oJq — uj 2 ) 3 and after their correction 
by including all couplings between adjacent oscillators (the coefficient iV/3 is 
substituted by N/2), we see that the results ( |10"1) and (|TT1) are different only in 
the coefficient Q\ + 2Q 2 cos9 which substitutes the coupling constant Qi in the 
eqs. (|7|) and (|8|). In the crystals belonging to the space groups of symmetry D% 
and .Df is cos 9 = —0.5 and the coefficient gives the value Qi — Q 2 . The couplings 
between even and odd oscillators have indeed the opposite efect on the ORD. 
Besides we include in our three coupled oscillators model all couplings between 
adjacent oscillators on the helices but only a half of couplings between even and 
odd oscillators. Nevertheless the influence of couplings with couplings constants 
Qi and Q 2 has the same form. If we include all couplings between even and odd 
oscillators the influence of these couplings will be, of course, much greater but 
we don't know the expression for this influence. 

Even and odd coupled oscillators on the helices are not contacted and we 
cannot use the method of the contacted compound oscillators from the section ^. 
If we want to solve this problem exactly we should solve the model of N coupled 
oscillators where N is the number of all single oscillators on one helix and in this 
model we should include all couplings between adjacent, even and odd oscillators. 
This model is not, of course, really solvable. 

From this reason we have proceeded by another way. We have solved further 
the models of four, five etc. coupled oscillators in which we have included the 
couplings between adjacent, even and odd oscillators and the results we have 
expressed only with the natural frequency u Q and only in the approximations of 
the oscillator strengths of the normal modes of vibrations and only for the crystals 
with the space groups of symmetry D\ and D|. Based on these results we have 
found the tendency in the expression that contains the coupling coefficient Q 2 and 
we have looked the limit of this expression with increased number of oscillators 
in the models. We show now the solving of the four coupled oscillators model. 
The solving of other models is similar but more complicated and we discuss only 
the results. 

3.1 The four coupled oscillator model of the optical rota- 
tory dispersion of crystals 

We extend the three coupled oscillators model by the fourth oscillator which 
lies on the helix in the plane z = 2d and his vibration direction is turned by the 
angle 29 around z axis with respect to the vibration direction of the oscillator 
in the plane z = 0. The coupling constant Q\ describes the coupling between 
adjacent oscillators and the constant Q 2 describes the couplings between the first 
and the third and between the second and the fourth oscillators. The motion 
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equations of the oscillators in the fields of the left and the right circularly polarized 
wave are 

h + wgn + Q x r 2 + Q 2 r 3 = F l { 



r 2 + uJqT 2 + Qin + Qtr 3 + Q 2 r 4 = F 2 ,r , 

h + ^o r 3 + Qir\ + Qi r 2 + Qiri = F 1 /, 



(12) 



r 4 + cjgr 4 + Q 2 r 2 + = F 1 / ', 



where r 1; r 2 , r 3 , r 4 are the displacements of oscillators from equilibrium, F denotes 
the projection of functioning forces to the motion direction of oscillators divided 
by the masses m of oscillators (electrons). In all terms in equation ([12]) and further 
in the sign ± hold the + sign for the left and the — sign for the right cicularly 
polarized wave. The electric field vector E has in our case the components 



(13) 



E l f = Eq cos (cut — ki >r z), 
E 1 ^ = ±E sin(ut — k\ r z) 

and then we can express the forces F^ r (r] = 1, 2, 3, 4) as 

F 1 / = ^ [(a cos9 + (3 sin 9) cos(ut + <j> ljf ) 
±(— a sin 9 + (3 cos 9) sin(u;t + 4>i,r)]i 
F 2 ' r = ^(a cos Lut±f3 sin ut) , 
F 1 / = *f*[(a£>s0 - /3sin0) cos(cut - <f> hr ) (14) 
±(a sin 9 + (3 cos 9) sm(ut — </>z, r )], 
F 1 / = ^ [(a cos 29 -p sin 29) cos(ut - 2(p hr ) 
±(asin 29 + (3 cos 29) sin{ut - 2<j> hr )\, 

where e is the electron charge and 0/ jr = k^ r = 27rra ' ,rd = ra ' ,r ^ d is a phase shift. 



We may express relations (|T2| ) in the normal coordinates which are 



qi = ^kW) {ri + Air2 + Air ' + Ul 
q 2 = -j^ TW) {r 1 + A 2 r 2 + A 2 r, + n\ 

Q^^^iri + A^-A^-n), 
g^^^in + A^-A^-n), 



(15) 



where 



Ai 
A 2 
A 3 
A, 



Qi + a/5Q?+8QiQ 2 - 


f-4Q| 


2(Qi+Q 2 ) 


Qi— V / 5Q?+8QiQ 2 - 


h4Q| 


2(Qi+Q 2 ) 






2(-Qi+Q 2 ) 


Qi + v /5Qf-8QiQ 2 - 


f4Q| 



2(-Qi+Q 2 ) 



(16) 
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As a result of coupling the natural frequency luq of each oscillator splits into four 
characteristic frequencies of the normal modes of vibrations 



u\ = uJq + 
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Qi-V 5 <2?+ 8< 3i £ 
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Qi-^5Qj-8Q 1 C 
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^3 = ^0 - 

, ,2 _ , ,2 

W 4 — ^ 2 

Now it is possible to rewrite the motion equations of each oscillator 

# r + w;gJ p = i2£, 77 = 1,2,3,4, (18) 

where i?^ are the forces in the normal coordinates. For them we get in general 
expression 

Rt = ^«) cosM + oj£), 77 = 1, 2, 3, 4 (19) 
where o l £ is a phase shift and for (a l £) we have derived 

KT) 2 = t 1 ± 3< ^V sin - 3 cos 9 + 4 cos 3 =p 120 1>P cos 2 OsinO 
— 2Ai (1 — cos 6* — 2 cos 2 ± 4>i, r sin $ =t 40z, r sin 9 cos 0) 
+A 2 (1 + cos 9 =F 0/,r sin 0)] 
(a^ r ) 2 = s£±|£ [1 ± 3</v sin - 3 cos 9 + 4 cos 3 9 t 120i, P cos 2 flsinfl 
—2A 2 (l — cos 9 — 2 cos 2 ± 0j )P sin^ ± 40; >r sin 9 cos 0) 
+A|(1 + cos 9 =F 0/, r sin 9)] 
(a^) 2 = [1 =|= 3</v sin + 3 cos - 4 cos 3 9 ± 120 1>P cos 2 9sin9 

+2A 3 ( 1 + cos 9 - 2 cos 2 9 T 4>i, r sin 6> ± 40^ sin 9 cos 9) 
+A 2 3 (1 - cos6» ± <f) l)r sm9)} 
(a^ r ) 2 = ^^±g [1 q= 30/, r sin + 3 cos 9 - 4 cos 3 9 ± 12^, r cos 2 0sin0 

+2A 4 (1 + cos 9 - 2 cos 2 9 =F 0/, r sin0 ± 40^ sin 6 1 cos 0) 
+A 2 (1 - cos0 ± ijr sin0)]. 



(20) 



Let's substitute eq. ([191) into eq. (18|), we get 



+ <4tf = «)^ cosM + < r ) (21) 



and the solution of these equations is 



, e_E n cos (tut + a l ' r ) 



The propagating light wave induces the dipole moments d l £ which are 

< = g„«)/ 9 „e (23) 
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or using 

J q e 2 E cosM + ^; 



9 " v 9 " ' m on 2 -on 2 v ' 

where f q are the oscillator strengths in the normal modes of vibrations. The 
mean polarizability per volume unit is 

N'd l f Nd l f 
x l f = *2 — — = *2 =— , (25) 

V E COS (cut + <Jql) 3E COs(ujt + (Tql) 

where N' is the number of compound oscilators in the volume unit. We see that 
in our case N' = N/3 where N is the number of isolated oscillators. 

For the refractive indices of the crystals we have the Drude-Sellmaier disper- 
sion relation 

4 AttNp 2 4 f 

n 2 r - 1 = AttYx 1 / = Yia'f) 2 ^—,. (26) 



From eq. (|2q) we are able to solve the relation 



2 „,_47rjvg A [«) 2 -K )U, 

3m ~, uj 2 — cu' 2 



nl -„:-- -. — Y ^^—^f^L (27) 

r,=l 



and using (|20| ) we can calculate for the expression in square brackets 

«) 2 - KJ 2 = ( " 2+ i+Af +<M sin 0(3-12 cos 2 9 

-2A l - 8 Ax cos 9 - A 2 ), 
KJ 2 - KJ 2 = ( " 2+/ ^ +<M sin 9(3 - 12 cos 2 

|2\ 



-2A 2 -8A 2 cos9-A 2 ), 
«) 2 - K 3 ) 2 = ( " 2+ f + ^ + ^ sinfl(-3 + 12 cos 2 9 
-2A 3 + 8A3Cos9 + A 2 ), 



(2i 



,1 v > _ i,,T y _ (a 2 +^)(^i+f) 

-2Aa + 8Aacos9 + A 2 ) 



l««r - l*«r = - -* + "g^ Binfl(-3 + 12 cos 2 9 



In eq. (|28|) we can write (<pi + <j> r ) = u)d ( ni + nr *> anc i n 2 _ n z = ( n ^ _ nr )( n ^ + 7^). 
Using the well known formula for the ORD p(u) = ^.( n i — n r ) we obtain 

. 2nNde 2 2 . 2 _ 2 . . . 

p(w) = ^ to 2 (a 2 + P 2 )sm9 



(3 - 12 cos 2 9 - 2At - 8A1 cos - A 2 )/, 



91 



(l + A 2 )(u 2 -u 2 ) 
(3 - 12 cos 2 9 - 2A 2 - 8A 2 cos - A 2 )/ 92 



(l+ADK-o; 2 
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(-3 + 12 cos 2 9 - 2A 3 + 8A 3 cos 9 + A 2 3 )f qs 
+ (1 + J®(ul - u?) 

(-3 + 12 cos 2 9 - 2A 4 + 8A 3 cos 9 + Aj)f q4 1 
+ (1 + A 2 )(^ 2 -^ 2 ) [ [J) 

Using eqs. ( |TTD we can rewrite the formula ( P5| ) to the form which contains 
the natural frequency of the oscillators only but we see that we would get the 
complicated expression. We will discuss the result ( |2"§| ) only in the approximations 
of the oscillator strengths and also only for the crystals with the space groups of 
symmmetry D| and D 3 for which 9 = 120 deg and therefore cos 9 = —1/2. 



3.2 Discussion 



We will express as in the sections [II |^ and [3] the result (29) in the linear 
harmonic oscillator approximation and in the Heitler - London approximation of 
the oscillator strengths. In all results we neglect all the small terms containing 
the expressions 1/(uJq—uj 2 ) of the order higher than second. These terms contain 
with the exception of the great value of the denominator also the higher orders 
of the small coupling constants Qi and Q2 in the numerators. 

At first we will use in the discussion of the result (^9[) the Chandrasekhar 
linear harmonic oscillator approximation where we can calculate with equation 
fqi — fq 2 = fqz = /?4 = /o- m this case we obtain the ORD formula of the 
Chandrasekhar type 



= — —-u^f- m 

The Heitler-London aproximation is in the case of four coupled oscillators 
given by the relation 

r f r f l 

(31) 



fqi fq2 fqz fqi fo 



UJ\ U0 2 UO3 UJ4 Uq 

and using eqs. (|TT| ) and taking again into account that Qi and Q2 are small 
quantities we can simplify expressions for the oscillator strengths to the form 



fqi 


= /o 


;i + 


/<?2 


= /o 


;i + 




= /o 


;i - 


fqi 


= /o 


1 - 



Qi + y/5Qj+8QiQ 2 +4Q 



Qi~y / 5Q^+8QiQ 2 +4Q 



4w 



2 1 



h—y/5Qf-8QiQ 2 +4Q| ' 



h + y/5Qi-&QiQ2+4:Q2 



Now we substitute these relations into eq. (|29| ) and we obtain 



27rNde 2 f 
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X 



3Qi - 4Q 2 2(3Q!CU 2 - 4Q 2 cu 2 - 2Q 2 + 2Q 1 Q 2 ) 



UJ 2 — UJ 2 (u) 2 — UJ 2 ) 2 



(33) 



and after neglecting the terms Q\ and Q\Q 2 and adding the terms in the square 
brackets we have the formula of the Agranovich type 

. 27rNde 2 f (a 2 + P 2 )sm9-(Q 1 -lQ 2 ) co 2 (co 2 + co 2 ) 

p{uj) = ^ ^M^y - (34) 

We see that applying the linear harmonic oscillator aproximation and the 
Heitler-London aproximation to the general ORD formula (p9|) of four coupled 
oscillator model we obtained for the crystals belonging to the space groups of 
symmetry D\ and the Chandrasekhar ([30]) and the Agranovich (|34|) formulae 
again. Both formulae contain the same expression Q\ — \Q 2 . Let's also note 



the results of solving three coupled oscillator model published in pC§ where these 
formulae contain the expression Q 1 — Q 2 . 

We obtain the similar results also in the five, six etc. oscillator models. In 
the linear harmonic oscillator and in the Heitler - London approximation we get 
again the Chandrasekhar and the Agranovich formulae which are different only 
in the constant containing the coupling constants Q\ and Q 2 . This constant has 
for the five oscillator model the form Qi — §Q 2 , in the six oscillator model the 
form Qi — |Q 2 etc. 

These results give the possibility to suppose the form of the ORD formulae 
for general N coupled oscillator model. We obtain the Chandrasekhar formula 
in the form 



AnNde 2 f (a 2 + (3 2 ) sm6 ■ (Q, - ^Q 2 ) w 2 

pH = ^ W -«?)*> (35) 

where A^ — 1 is the number of coupling between neighbouring oscillators and N — 2 
the number of coupling between odd and even oscillators in the model. It means 
that for the practical case N — > 00 we have 

. A7cNde 2 (a 2 + p 2 )sm9-(Q 1 -2Q 2 ) lu 2 , . 

P{UJ) = -7-2 =rj. (36) 

mC 2 [UJq — UJ 2 ) 2 

The same results we can write for the Agranovich formula: 

27rNde 2 f (a 2 + {3 2 )sm6-(Q 1 - 2 -^Q 2 ) to 2 (u 2 + u 2 ) 

p(uj) = ^ ugK 2 -" 2 ) 2 m 

and for N — > 00 

. 27rNde 2 f (a 2 + [3 2 )sm9-(Q 1 -2Q 2 ) u 2 (oo 2 + u 2 ) 

p(w) = ^ u; >g ( } 

We see that the expression Qi — 2Q 2 is the common limit of the coefficients 
containing the coupling constants in the final ORD formulae in the case N — > 00. 
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4 Conclusions 



In this paper we have proved that with the Chandrasekhar ORD two coupled 
oscillator model we can include all couplings between adjacent oscillators on the 
helix that is given by the crystal structure. We have solved that the influence 
of the couplings between adjacent oscillators in the cases when the second single 
oscillator in the first compound oscillator is the first oscillator in the second 
compound oscillator we can hold as aditive (the mistake in this conclusion is 
neglectable). For the including of all couplings between adjacent oscillators is 
sufficient to take in the Chandrasekhar model that the number of single oscillators 
in the volume unit is the same as the number of the couplings between adjacent 
oscillators. The same results were proved for the models of three and four coupled 
oscillators - the parts of the ORD results that contain the coupling constant 
between adjacent oscillators Q\ are identical in all models. 

From the crystal structure it follows that we cannot neglect the couplings 
between even and odd oscillators on the helix. We can assume for the real crystal 
that the value of the coupling constant between even and odd oscillators Q2 is 
comparable with the value of the coupling constant between adjacent oscillators 
Qi. Further we have proved that in the crystals with the space groups of symme- 
try D\ and the couplings between even and odd oscillators have the opposite 
effect on the ORD than the couplings between adjacent oscillators. 

The first coupled oscillator model in which we can include the coupling be- 
tween even or odd oscillators is the three oscillator model. But using the previous 
conclusion this model can include all couplings between adjacent oscillators and 
only a half of couplings between even and odd oscillators on the helix. Besides we 
cannot assume these couplings as aditive because they function between single 
oscillators inside of the compound oscillators - for example in the three oscillator 
model the second oscillator in the compound oscillator is in the fact coupled also 
with the second oscillator in the adjacent compound oscillator (becouse in the 
used model the last single oscillator in any compound oscillator is the first single 
oscillator in the adjacent compound oscillator). 

From this reason we have further solved the models of four, five etc. coupled 
oscillators. In these models the number of the included couplings between even 
and odd oscillators increases and the results can be generalized. By this way 
we have proved that the couplings between even and odd oscillators also have 
the aditive effects. The ORD results of the three oscillators model after the 
approximations of the oscillator strengths of the normal modes of vibrations 
contain the coefficient Q± — Q2 (for the crystals with space groups of symmetry 
D\ and D®). In the limit of N coupled oscillators for N — > 00 it was derived the 
coefficient Q\ — 1Q 2 . This also means that the relative influence of the couplings 
between even and odd oscillators on the ORD is twofold in comparison with the 
couplings between adjacent oscillators. 

We can formally obtain the generalized ORD results of the N coupled oscil- 
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lators model after the approximations of the oscillator strengths using the three 
oscillators model in which we denote the value of coupling between the first and 
the third oscillators as 2Q2. We can say that this conclusion is practically ac- 
ceptable only after using the results of the three oscillator model in the numerical 
approximations of the experimental ORD data for example of a-quartz or tellur. 
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